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1.  Magnetic  Field  and  Particle  Flux  Specification  in  the  MSM/MSFM 


1.1.  Original  Mapping  Method 

T\ytHilmer  and  Voigt  [1995]  magnetic  field  model  provides  the  Magnetospheric  Specification 
Model  (MSM)  [Freeman  et  al,  1993]  and  Magnetospheric  Specification  and  Forecast  Model 
(MSFM)  [Bales  et  al,  1993]  with  equatorial  magnetic  field  values,  magnetic  flux  tube  volumes, 
and  a  method  for  establishing  one-to-one  mappings  between  an  ionospheric  electric  potential  grid 
and  a  simulation  grid  contained  within  the  GSM  equatorial  plane.  Note  that  the  MSM  and  MSFM 
utilize  identical  magnetic  field  treatments  in  all  respects  so  all  statements  and  conclusions  herein 
apply  equally  to  both.  The  MSM  simulation  grid  is  always  aligned  with  the  GSM  equatorial  plane 
owing  to  the  zero  dipole  tilt  approximation  which  was  adopted  to  take  advantage  of  North-South 
magnetic  field  symmetries  and  reduce  the  amount  of  pre-calculated  magnetic  field  configuration 
information  tabulated  and  stored  off-line.  Allowing  the  dipole  tilt  angle,  y/,  to  vary  would 
introduce  a  fourth  parameter  space  variable  (the  others  being  Dst,  the  magnetopause  standoff 
distance,  and  the  midnight  equatorward  auroral  boundary)  and  would  increase  the  number  of 
magnetospheric  configurations  by  at  least  a  factor  of  five.  The  MSM  algorithm  utilizes  stored 
magnetic  field  quantities  because  the  on-line  tracing  of  magnetic  field  lines  is  too  computer 
intensive  to  be  practical  in  an  operational  setting.  These  magnetic  field  quantities  are  used  by  the 
MSM  to  map  the  electric  field  (magnetic  field  lines  are  assumed  to  be  equipotentials)  to  the 
equatorial  plane  and  calculate  adiabatic  particle  drift  paths.  By  performing  drift  calculations  on  a 
variety  of  particle  species  over  appropriate  energy  ranges,  one  of  the  MSM  final  outputs  is  the 
specification  of  particle  fluxes  in  the  magnetic  equatorial  plane. 

A  special  procedure  developed  by  Hilmer  et  al  [1993],  commonly  referred  to  by  its 
FORTRAN  subroutine  name  MAP3D,  utilizes  additional  tabulated  magnetic  field  information  to 
obtain  MSM  particle  flux  and  magnetic  field  information  outside  the  GSM  equatorial  plane.  In 
the  MAP3D  procedure,  a  3-D  grid  was  defined  and  magnetic  field  vector  and  equatorial  mapping 
information  stored  for  each  grid  point  in  each  of  932  magnetospheric  configurations  spanning  the 
3-D  parameter  space.  To  represent  actual  geophysical  conditions,  the  MAP3D  algorithm 
interpolates  magnetic  field  information  from  the  eight  closest  configurations  in  the  parameter 
space  to  form  a  unique  configuration  for  each  given  time.  Arbitrary  points  in  space  are  then 
mapped  to  the  MSM  simulation  region  using  this  tailored  field  configuration  and  the 
corresponding  equatorial  fluxes  are  scaled  (using  the  FLUX3D  algorithm)  to  assign  appropriate 
values  to  off-equatorial  locations.  As  with  the  MSM  particle  drift  algorithm  mentioned  above, 
this  procedure  takes  advantage  of  the  symmetry  inherent  in  the  magnetic  field  model  when  the  tilt 
of  the  Earth's  dipole  field  is  fixed  equal  to  zero  degrees. 
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1.2.  Mapping  Problems 

With  the  established  procedure  of  providing  location  input  in  GSM  coordinates  regardless  of 
the  actual  dipole  tilt  angle,  the  MAP3D  algorithm  has  been  used  to  extract  reasonable  particle  flux 
values  from  the  2-D  MSM  simulation  output  under  a  variety  of  geophysical  conditions  and  along 
very  different  orbits  including  the  high  inclination  CREIES  orbit.  Recently,  however,  it  was 
revealed  that  some  unphysical  magnetic  field  mappings  were  being  produced  by  MAP3D  with  the 
implication  that  the  MSM  particle  flux  specifications  would  be  adversely  affected  (Prochaska, 
Hughes  STX  Corporation  -  23  April  1996  facsimile  to  Freeman  of  Rice  University).  Prochaska 
noticed  on  several  different  occasions  that  magnetic  field  lines  from  geosynchronous  satellites  on 
the  nightside  of  the  magnetosphere  mapped  to  equatorial  locations  far  tailward  of  geosynchronous 
altitudes,  i.e.,  beyond  “  “20  Re  (Earth  radii).  This  observation  is  most  certainly  indicative 
of  a  mapping  problem  as  all  dipole  tilt-dependent  magnetic  field  models  we  know  of,  including  the 
Hilmer  and  Voigt  [1995]  model  on  which  MAP3D  is  based,  map  geosynchronous  orbit  locations 
to  neutral  sheet  locations  within  just  a  few  Re  of  the  satellite's  position.  Several  causal  factors 
were  suggested  including  the  presence  of  extremely  stretched  magnetic  field  lines  in  the  midnight 
sector  of  the  magnetotail  and  large  satellite  magnetic  latitude  deviations  (>  1 1  degrees)  as  a 
function  of  geographic  longitude.  While  both  of  these  certainly  contribute  to  magnetic  field 
mappings  beyond  geosynchronous  altitude,  the  overwhelming  contributing  factor  involves  the  y/~ 
0  approximation  adopted  within  the  MAP3D  algorithm.  Geosynchronous  mapping  examples 
illustrates  this  point. 

First,  we  ran  the  MSM  for  a  short  period  with  .^  =  3  in  order  to  use  the  MAP3D  algorithm 
at  UT=  0500  on  day  354  of  1995  when  the  actual  dipole  tilt  angle  was  The 

coordinates  (Jl'(,^,yp5^^,ZQ5j^)  =  (-6.065Re, -0.123 Re,-2.626Re)  represent  a  point  along  the 

geosynchronous  orbit  at  that  time.  With  the  y/=0  approximation,  the  eight  magnetic  field 
configurations  interpolated  by  MAP3D  scattered  mappings  from  =-9  to  -37  Re  such 
that  the  resultant  mapping  was  located  well  down  the  magnetotail  at  X^^  =  -27  Re.  Second,  if 
this  same  mapping  is  done  using  the  Kp  =  3  version  of  the  Hilmer  and  Voigt  [1995]  model  with 
the  dipole  tilt  angle  also  held  at  y/=0,  this  geosynchronous  point  maps  to  X(.^  =  -33  Re.  In 
both  ^  =  0  cases,  the  starting  point  is  located  on  moderately  stretched  tail-like  field  lines  more 
than  2  Re  below  the  magnetic  neutral  sheet  (and  the  MSM  simulation  region).  The  neutral  sheet 
is  the  surface  defined  by  the  local  magnetic  field  strength  minimum  center  in  the  plasma  sheet  and 
coincides  with  the  GSM  equatorial  plane  when  y/=  0.  In  fact,  locations  near  the  center  of  the 
magnetotail  where  jZ^^I  >  1.5  to  2.0  Re  and  A’Q^=-6.6Re  often  map  well  tailward  of 

geosynchronous  orbit  when  ICp  is  moderate  to  high  and  Finally,  if  the  actual  value  of 

-34°  is  used  in  the  Hilmer  and  Voigt  [1995]  model  the  magnetic  neutral  sheet  is  displaced 
below  the  GSM  equatorial  plane  and  our  geosynchronous  point  maps  more  realistically  to  the 
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neutral  sheet  inside  of  =  -7  Re.  This  close  mapping  results  directly  from  that  fact  that  our 
point  is  closer  to  the  neutral  sheet  which  has,  in  this  latter  case,  been  properly  displaced  out  of  the 
GSM  equatorial  plane  in  response  to  a  non-zero  dipole  tilt  angle. 

We  conclude  that  the  implementation  of  the  MAP3D  algorithm  must  be  modified  to 
incorporate  the  effects  of  dipole  tilt  angle  variations  on  the  magnetic  field. 

1.3.  Proposed  Solution 

With  ?/  =  0,  the  GSM  equatorial  plane  is  coplanar  with  the  MSM  simulation  region, 
extending  approximately  10  Re  sunward  and  20  Re  tailward  of  the  Earth,  and  acts  as  a  symmetry 
plane  for  the  MAP3D  field  mapping  algorithm.  This  approximation  causes  the  plasma  sheet  and 
its  imbedded  magnetic  neutral  sheet  to  be  artificially  aligned  with  the  GSM  equatorial  plane.  To 
be  more  physical,  the  flat  MSM  simulation  region  should  flex  and  follow  the  motion  of  the  plasma 
and  the  magnetic  neutral  sheets  as  they  are  displaced  from  the  GSM  equatorial  plane  in  response 
to  phenomena  such  as  dipole  tilt  angle  changes,  solar  wind  pressure  variations,  and  the  growth 
and  decay  of  the  ring  and  cross-tail  currents.  Accordingly,  the  MAP3D  algorithm  should 
incorporate  the  corresponding  magnetic  field  configuration  changes.  As  both  of  these  alteration 
would  require  much  additional  research,  testing,  and  computer  storage  with  no  guarantee  of 
improving  MSM  particle  flux  specifications  significantly,  we  propose  a  relatively  simple  solution 
in  the  form  of  a  coordinate  transformation  that  will  allow  the  MSM  and  MAP3D  algorithms  to 
remain  imchanged. 

As  mentioned  above,  while  moving  away  from  the  neutral  sheet  of  the  magnetotail  along  the 
^GSki  direction  we  cross  field  lines  that  map  increasing  farther  tailward.  Even  within  the  MSM  and 
MAP3D  framework,  this  fact  means  that  we  are  encountering  particles  with  gradually  more 
distant  neutral  sheet  crossing  points  and  the  particle  fluxes  extracted  from  MSM  simulations  using 
MAP3D  are  related  to  the  flux  levels  at  those  same  points.  We  are  currently  using  the  GSM 
coordinates  of  satellites  to  derive  MSM  flux  levels  using  MAP3D  regardless  of  the  actual  neutral 
sheet  geometry.  We  need  to  modify  the  input  position  information  provided  to  MAP3D  to 
account  for  realistic  motions  of  the  neutral  sheet. 

To  improve  the  magnetic  field  mappings  used  to  derive  particle  flux  specifications  from  both 
MSM  and  MSFM  simulations,  we  propose  providing  the  MAP3D  algorithm  with  input 
coordinates  based  on  a  new  magnetospheric  neutral  sheet-oriented  coordinate  system  dependent 
on  both  the  dipole  tilt  angle,  yr,  and  geomagnetic  activity,  Kp. 
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2.  A  Neutral  Sheet-Oriented  Coordinate  System 


2.1.  Characteristics  of  the  Neutral  Sheet 

For  the  purposes  of  this  work,  we  identify  the  magnetospheric  neutral  sheet  simply  as  the 
minimum  magnetic  field  |B|  surface  found  at  low  latitudes  about  the  Earth  at  all  local  times. 
When  the  dipole  tilt  angle  y/\&  zero  the  magnetic  field  configuration  becomes  symmetric  about  the 
^GSM  -  ®  plane  and  the  neutral  sheet  is  coincident  with  the  GSM  equatorial  plane.  When  y/^0, 
this  alignment  does  not  occur.  For  positive  (negative)  y/,  the  neutral  sheet  is  shifted  above 
(below)  the  GSM  equatorial  plane  at  midnight  as  represented  in  Figure  1  (with  \}/=  i/).  The 

neutral  sheet  bends  gradually  away  fi'om  the  SM  equator  to  become  aligned  with  the  axis  at 
Zqsm  =  -HD  sin  y/,  where  HD  is  known  as  the  "hinging"  distance,  and  there  is  no  tendency  for  it  to 
return  to  the  GSM  equatorial  plane  at  large  distances  {Fairfield,  1980].  As  we  move  toward  the 
flanks  of  the  mid  and  far  magnetotail  the  neutral  sheet  approaches  the  GSM  equatorial  plane  and 
crosses  it  [e.g.,  Fairfield,  1980;  Dandouras,  1988].  Away  fi’om  local  midnight,  the  near-Earth 
neutral  sheet's  deflection  fi’om  the  SM  equatorial  plane  lessens  as  it  joins  the  SM  equatorial  plan  at 
dawn  and  dusk  [Lopez,  1990].  Finally,  on  the  dayside  we  assume  for  simplicity  that  the  neutral 
sheet  remains  fixed  in  the  SM  equatorial  plane  which  contains  the  local  |B|  minimum  of  the  dipole 
field.  As  these  various  neutral  sheet  regions  must  connect  smoothly  to  form  a  single  surface,  the 
neutral  sheet  must  be  continuously  flexing  and  warping  in  response  to  the  diurnal  dipole  tilt  angle 
oscillations.  Additionally,  there  is  a  tendency  for  the  nightside  neutral  sheet  to  be  displaced  more 
from  the  GSM  equatorial  plane  as  geomagnetic  activity  decreases  for  a  given  value  of  y/  (see 
details  in  Section  2.3),  a  feature  that  is  as  least  partially  attributable  to  the  existence  of  weaker 
cross-tail  currents  during  these  times.  Without  detailed  quantitative  information  about  how  to 
model  the  neutral  sheet's  response  to  other  more  specific  physical  phenomena,  e.g.,  solar  v^nd 
pressure  changes,  we  will  concentrate  on  including  an  activity  dependence  via  the  geomagnetic 
index  in  hopes  that  it  will  aide  us  in  specifying  neutral  sheet  behavior. 

In  the  following,  we  derive  a  single  analytic  coordinate  system  {X',Y',Z')  dependent  on  y/ 

and  Kp  such  that  the  Z'  =  0  surface  coincides  approximately  ’with  the  neutral  sheet.  The  GSM 
coordinates  of  satellite  locations  used  previously  to  determine  flux  specification  should  be 
translated  into  this  new  coordinate  system  before  being  supplied  to  the  MAP3D  algorithm. 

2.2.  Dependence  on  Dipole  Tilt  Angle 

The  Coordinates  X'  and  Z.  The  geometry  used  to  organize  the  new  coordinate  system  in  a 
given  GSM  X-Z  plane  is  shown  in  Figure  1 .  The  SM  coordinate  system  differs  from  GSM  by  a 
rotation  equal  to  the  dipole  tilt  angle  about  the  common  =  Y^  axis,  namely 
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Neutral  Sheet-Oriented  Coordinate  System 


Schematic  of  neutral  sheet-oriented  coordinate  system  regions  and  boundaries  in  a  GSM  X-Z 
plane  (Sun  to  the  left).  Contours  of  constant  Z'  are  equidistant  from  the  neutral  sheet  (at  Z'  =  0) 
while  JT  increases  along  them  to  the  left  and  T  increases  out  of  the  page.  The  dipole  tilt  angle  is 
^/and  the  neutral  sheet  "hinging"  distance,  HD,  varies  with  Kp.  The  angle  if/,  equal  to  y/  at  T^w 
=  0,  decreases  in  magnitude  toward  the  flanks  and  |</|  <  |^^^.  The  point  (XC,  ZC)  acts  as  the 
center  of  curvature  for  the  curved  coordinate  space  in  Region  R2.  Restrictive  conditions  include 
I  ^  <  40  degrees  and  IZ^^J  <Z0=  17  Re. 


FIGURE  1 
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^SM  ^GSM  '^GSM  W 


iw = (i) 

^SM  =  -^GwSin 

We  define  the  new  SM'  coordinate  system  similarly  except  that  its  rotation  angle  y/  varies  with 
^Gw  ^  modified  version  of  Eq.  (10a)  of  Dandouras  [1988],  which  incorporates  a  nominal 
"hinging"  distance  of  10.5  Re,  is  used  to  get  the  position  of  the  neutral  sheet,  Z^^^y  across  the 
magnetotail  for  <  XC  and  illustrates  the  relationship  between  the  angle  y/  and  namely 

HD  (  Y 

^ ^  <  13.86Re 

(2) 

^ GSM. NS  I^GSw]  —  13.86Re 


The  angle  y/  is  largest  and  equal  to  y/^  at  the  center  of  the  magnetotail  and  decreases  toward  the 
flanks  to  y/  =  0  for  |  >  13.86  Re.  The  Z'  =  0  contour  is  the  line  labeled  as  the  neutral  sheet 

while  a  sample  negative  Z  contour  is  indicated  by  the  line  running  parallel  to  it.  For  \Zq^  \<Z0, 
each  plane  is  divided  into  Regions  RO,  Rl,  R2,  and  i?5,  separated  by  the  three  lines  emanating 
fi'om  the  center  of  curvature  point  (XC,  ZC).  The  constant  ZO  is  selected  to  be  just  earthward  of 
the  smallest  \ZC\  encountered.  The  Region  RO-Rl  boundary  is  parallel  to  the  Z^  axis,  the  Region 
R1-R2  boundary  passes  through  the  Y^^  axis,  and  the  Region  R2-R3  boundary  is  parallel  to  the 
Zq^  axis.  As  yr  changes  sign  fi'om  positive  to  negative  the  neutral  sheet  in  the  magnetotail  shifts 
fi'om  above  to  below  the  GSM  equatorial  plane  by  an  amount  related  to  the  "hinging"  distance 
HD  and  the  angle  y/.  Meanwhile,  the  neutral  sheet  on  the  dayside  remains  aligned  with  the  SM 
equatorial  plane  and  is  displaced  in  the  opposite  sense.  We  now  define  of  the  neutral  sheet- 
oriented  coordinates  X,  Y',  and  Z  in  each  of  the  four  regions. 

Region  RO.  In  Region  RO,  the  coordinates  Z  and  Z^-^  are  equivalent  while  X  is  aligned  with 
and  varies  along  the  X^^  axis  as  a  function  of  the  difference  between  y/  and  y/,  namely 
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such  that  the  SM*  and  SM  coordinates  are  equivalent  at  =  0  where  we  define  \f/=  y/' . 

Region  Rl:  In  Region  RI,  the  JC  coordinate  is  defined  as  in  Region  RO  but  the  Z*  coordinate 
is  modified  to  account  for  this  region's  wedge  shape  and  written  as 


Z'  =  - 


cos(5</-  -  y/') 


and  Z'  ■■ 


'SM 


cos{y/-y/'  -0) 


(4) 


where  9  is  the  angle  between  the  Region  RI-R2  boundary  and  the  point  of  interest  in  Region  RI 
measured  about  the  point  {XC,  ZQ.  This  provide  a  smooth  transition  in  Z  across  the  Region  RO- 
R1  boundary  where 

Region  R2:  Circular  arcs  measured  about  (XC,  ZQ  are  used  to  define  constant  Z  contours 
such  that  the  zero  contour  passes  through  the  GSM  origin  at  the  earthward  boundary  of  Region 
R2.  The  Region  R2-R3  boundary  is  the  vertical  line  where  X^^  =  XC .  If  Rq  is  the  distance 
fi-om  (XC,  ZQ  to  some  point  of  interest  Region  R2,  then  a  point  (^,2,212)  on 

the  Region  RI-R2  boundary  with  the  same  Z'  value  has  the  coordinates 


X„  =  XC  +  sin  y/ 


=  ZC  +  -^Rc  cos  y/' 


(5) 


Our  Region  R2  neutral  sheet-oriented  coordinates  can  then  be  written  as 


X'  =  -Rc  cos' 


r 


1.0-0.5 


R. 


[Rc-\HDsiny/'-ZC^ 


(6) 


where  X  is  the  distance,  measured  along  a  circular  arc  of  radius  Rq,  between  the  Region  R1-R2 
boundary  point  {X,^,Z,^)  and  any  Region  R2  point  (Xoa^,ZGaf).  We  see  that  X  is  always 

negative  in  Region  R2  while  Z  is  positive  above  and  negative  below  the  neutral  sheet. 
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Region  R3:  The  Region  R3  system  becomes  rectilinear  and  its  coordinates  are  written  as 


X' = -|(z„w  -  zcy  I  - 

^'  =  ^GSM  -HDsmy/' 


(7) 


The  coordinate  X,  measured  tailward  from  the  Region  R1-R2  boundary,  incorporates  an 
appropriate  arc  length  from  Region  R2  and  a  straight  portion  tailward  of  X^^  =  XC  while  X  is 
measured  along  the  direction  relative  to  the  neutral  sheet  position. 

The  Coordinate  ¥':  Our  new  coordinate  T  depends  on  the  value  of  the  "hinging" 
distanced),  and  the  dipole  tilt  angle  yr.  We  integrate  Eq.  (2)  to  get  the  length  along  the  curve  it 
defines  in  the  GSM  Y-Z  place  between  =  0  and  the  point  of  interest  and  get 


Y'  =  +  (l  +  ^2  y J|  I  <  1 3 .86  Re 

where  =  (l  •439 1 534  x  1 0”^  )sin  V 


(8) 


When  outside  the  indicated  range  we  get  Y  by  evaluating  (8)  at  Y^^m  =  ±13.86  Re,  as 
appropriate,  and  adding  a  length  equal  to  ^Y^^m  \  -13.86)  Re  owing  to  the  linear  flank  portion. 

This  flattening  of  the  neutral  sheet  in  the  outer  magnetotail  flanks  will  not  significantly  affect 
magnetic  field  mappings  as  the  observed  neutral  sheet  deviates  only  slightly  (in  the  opposite  sense 
from  the  midnight  deflection)  from  the  GSM  equatorial  plane  in  that  region. 

2.3.  Dependence  on  Magnetic  Activity 

The  coordinate  geometry  described  above  requires  knowledge  of  the  neutral  sheet  "hinging" 
distance  HD  which  is  the  geocentric  distance  measured  along  the  negative  SM  x-axis 
characterizing  the  neutral  sheet's  displacement  from  the  =  0,  namely  Az  =  HD  sin  yr.  As 

noted  above,  most  neutral  sheet  models  have  adopted  the  use  of  a  single  average  value  for  HD 
and  thus  are  incapable  of  representing  shifts  related  to  variations  in  cross-tail  current  strength  and 
distribution.  In  magnetic  field  models,  for  example,  it  can  easily  be  shown  that  decreasing  the 
cross-tail  current  density  increases  the  relative  influence  of  the  tilting  main  field  so  the  neutral 
sheet  is  displaced  farther  from  the  equatorial  plane.  As  reviewed  briefly  by  Fairfield  [1980], 
observations  clearly  indicate  the  tendency  for  the  neutral  sheet  to  be  closer  to  the  GSM  equatorial 
during  magnetically  active  times  (when  the  cross-tail  current  density  is  larger  on  average).  This 
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implies  that  our  coordinate  system  geometry  would  benefit  if  HD  were  dependent  on  magnetic 
activity,  i.e.,  if  we  had  HD  =  HD{Kp). 

As  a  guide,  we  will  use  the  relationship  Lopez  [1990]  determined  using  AMPTE/CCE 
magnetic  field  data  firom  the  nightside  region  between  R  =  5  and  8.8  Re,  namely 

A4L^7’  =  -(0.14^;?  +  0.69)[cos<I»f^(0.065i?°*  -0.16^  (9) 


where  MLAT  is  the  position  in  degrees  of  the  neutral  sheet  relative  to  the  SM  equator,  Kp  is  the 
magnetic  activity  index,  3>  is  the  magnetic  local  time  in  degrees  (O  =  0  at  midnight),  R  is 
geocentric  radius  in  Re,  and  t^is  the  dipole  tilt  angle  in  degrees.  Eq.  (9)  describes  a  neutral  sheet 
that  diverges  increasingly  fi'om  the  SM  equatorial  plane  as  each  of  the  variables  R,  and  y/ 
increase  in  magnitude  while  the  neutral  sheet  converges  with  the  equatorial  plane  at  both  dawn 
and  dusk.  This  latter  feature  supports  our  selection  of  the  SM  equatorial  plane  for  the  neutral 
sheet  position  on  the  dayside. 

An  expression  for  HD  as  a  fiinction  o^Kp  was  derived  by  comparing  the  neutral  sheet  position 
estimates  given  by  Eq.  (9)  with  those  of  our  new  coordinate  system.  Using  a  large  dipole  tilt 
angles  of  yr=  +34°,  the  neutral  sheet  position  given  by  (9)  was  determined  for  the  full  range  of 
Kp  values  (Kp  =  0.0  to  9.0  by  0.3333)  and  R  values  {R  =  5.0  to  8.8  Re  by  0.1  Re)  at  midnight. 
The  equivalent  coordinates  in  our  new  neutral  sheet  oriented  coordinate  system  are  found  and  the 
error  between  the  two  methods  is  represented  by  the  magnitude  of  our  new  Z'  coordinate  as 
ideally  the  Z'=0  surface  coincides  with  the  neutral  sheet.  The  linear  function  minimizing  the  error 
is  given  by  the  upper  expression  of  the  following  equation  used  for  the  entire  Kp  range. 


HD  =  {-QM)Kp +  1236661  for  6.0  <Kp<  6.0 
HD  =  8.20667  for  Kp  >  6.0 


(10) 


where  the  "hinging"  distance  HD  is  in  Re  and  ranges  fi'om  approximately  13.4  to  5.6  Re  as  Kp 
increases  from  values  0.0  to  9.0.  The  maximum  error  (given  by  \Z’\),  occurring  for  the  extreme 

Kp  values  of  0.0  and  9.0,  is  0.133  Re  while  the  average  error  approaches  zero  (0.0004  Re)  for  a 
value  ofKp  very  close  to  3.33  (or  3+)  whenHD  equals  its  nominal  value  of  10.5  Re.  The  average 
overall  error  is  0.086  Re  for  this  y/=  ±34°  case,  is  half  as  large  when  ±19°,  and  vanishes  as 
expected  when  0°.  Because  GSM  input  is  restricted  to  \Zq^  |  <  ZO  and  ZO  depends  on  the 

value  of  HD  (relationship  not  shown  here)  we  must  restrict  use  of  the  above  linear  relationship  to 
Kp  values  less  than  6.0  in  order  to  allow  {Z^^J  input  values  of  up  to  ZO  =  17  Re.  This 
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modification  introduces  some  additional  error  only  in  the  high  Kp  range  but  increases  the  overall 
average  error  fi’om  0.048  to  0.068  Re  at  midnight  for  ±19°  case.  Faiifield  [1980]  notes 

that  there  are  no  large  changes  in  erroneous  neutral  sheet  prediction  when  our  range  of  HD 
values,  8.2  to  13.4  Re,  is  used.  We  note  that  while  the  error  has  been  minimized  at  midnight  and 
is,  by  design,  zero  at  dawn  and  dusk,  the  average  error  between  our  neutral  sheet  model  and  Eq. 
(9)  peaks  at  0.13  Re  near  MLT  values  of  5  and  19  for  the  ±19°  case.  These  differences  are 
related  to  the  dependence  of  the  ellipse  used  to  describe  the  far-tail  neutral  sheet  position. 

3.  Examples 

Figures  2  to  39  show  neutral  sheet-oriented  coordinate  {X,  Y,  Z)  grids  plotted  in  various  X-Z 
and  Y-Z  GSM  planes  for  conditions  defined  by  permutations  of  a  selection  of  y/  and  Kp  values 
(i.e.,  y/=  5,  20,  and  35  degrees  and  Kp  =  0,3,  6).  Contours  of  constant  X,  F,  and  Z  assume  the 
same  approximate  relation  to  each  other  as  the  straight  contours  of  constant  X,  Y,  Z  do  in  the 
GSM  coordinate  system.  In  each  plot,  the  Earth  is  represented  by  a  cross  at  the  projected 
location  of  the  GSM  origin.  For  y/=0  degrees,  the  neutral  sheet-oriented  coordinates  becomes 
equivalent  to  the  GSM  coordinates. 

In  the  Xq^-Zq^  grid  views  (with  Sun  to  the  right),  the  vertices  formed  by  intersecting 
contour  lines  remain  equidistant  as  the  grid  flexes  with  changing  dipole  tilt  angle,  y/.  In  Region 
R2,  Z'  contours  form  circular  arcs  centered  at  the  point  (XC,  ZC)  and  are  tangent  to  the  X'^ 
direction  at  the  Region  R1-R2  boundary  and  tangent  to  the  X^^  direction  at  the  Region  R2-R3 
boundary.  See  Figure  1  (with  Sun  to  the  left)  for  Region  and  boundary  definitions.  The  defined 
neutral  sheet  represents  the  Z'  =  0  surface  such  that  the  constant  Z  contour  highlighted  in  Figure 
1  has  a  value  Z'  <  0.  The  X  coordinate  values  are  positive  sunward  of  the  X'^  =0  plane.  In 

Region  R3,  the  neutral  sheet  is  displaced  in  the  Z^^m  direction  by  an  amount  proportional  to  the 
neutral  sheet  "hinging"  distance  HD,  which  decreases  as  Kp  increases.  For  positive  (negative)  tilt 
angles  y,  the  displacement  is  above  (below)  the  GSM  equatorial  plane. 

In  the  Y^^-Zfj^  views,  grid  vertices  formed  by  intersecting  contour  lines  are  equidistant  only 
for  Xq^  <  XC  as  contours  of  constant  Y'  are  defined  to  be  straight  lines.  Down  the  magnetotail 
center,  e.g.,  at  Xq^  =  -21  Re,  the  Z'  contours  become  elliptic  as  the  neutral  sheet  is  displaced 
above  (below)  the  GSM  equatorial  plane  for  positive  (negative)  y/  values.  Moving  toward  the 
magnetotail  flanks,  the  neutral  sheet  approaches  and  becomes  aligned  with  the  GSM  equatorial 
plane.  As  we  move  sunward,  i.e.,  from  X^^  =  -7  Re  to  +7  Re,  the  grid  patterns  appear  more 
irregular  as  the  selected  viewing  plane  intersects  multiple  Regions.  On  the  dayside,  Z'  =  0 
contours  are  straight  lines  in  the  SM  equatorial  plane. 
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4.  Comments 


The  new  magnetospheric  coordinate  system  presented  here  responds  to  changes  in  the  Earth's 
dipole  tilt  angle  and  the  Kp  index.  By  tracking  the  motion  of  the  magnetic  neutral  sheet  as  a 
function  of  time  and  geomagnetic  activity,  it  can  be  used  to  estimate  the  relative  position  of 
arbitrary  points  with  respect  to  this  magnetic  field  minimum  surface.  Accurate  knowledge  of  the 
neutral  sheet's  position  is  important  because  it  greatly  influences  the  shape  and  tailward  extent  of 
magnetic  field  lines  connecting  the  ionosphere  with  the  plasma  sheet.  The  effectiveness  of  the 
MAP3D  procedure  used  to  determine  energetic  particle  fluxes  for  arbitrary  locations  outside  the 
MSM  and  MSFM  simulation  planes  is  compromised  because  it  does  not  account  for  variations  in 
neutral  sheet  position,  i.e.,  the  geomagnetic  field  configurations  it  uses  fix  the  tilt  angle  equal  to 
zero.  It  is  suggested  that  improved  particle  flux  specifications,  owing  to  more  realistic  magnetic 
field  mapping  assignments  to  the  simulation  planes,  will  result  if  all  future  spatial  inputs  to  the 
MAP3D  algorithm  are  expressed  in  this  new  coordinate  system  rather  than  in  GSM  coordinates. 
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FIGURE  10 
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